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Abstract. We prove a law of large numbers and a central limit theorem for 
a tagged particle in a symmetric simple exclusion process in Z with variable 
diffusion coefficient. The scaling limits are obtained from a similar result for 
the current through —1/2 for a zero-range process with bond disorder. For the 
CLT, we prove convergence to a fractional Brownian motion of Hurst exponent 
1/4. 



1. Introduction 

A classical problem in statistical mechanics consists in obtaining the scaling limit 
of a tagged (sometimes called tracer) particle in a system of interacting particles. 
In general, we expect the limiting process to be a Brownian motion for the system 
in equilibrium |15| . or a timc-inhomogencous diffusion in general |13| . However, 
in the special case of a one-dimensional, nearest-neighbors system (the so-called 
single- file diffusion) , a sub-diffusive behavior is expected [9] . In the case of a sym- 
metric, simple exclusion process, the scaling limit of a tagged particle is given by 
■\/(l — ip)/ip)Zt, where Zt is a fractional Brownian motion of Hurst exponent 1/4 
and ip is the density of particles [1] , [5] . 

In recent years, the evolution of a random walk in random environment have 
received a lot of attention. Therefore, to investigate the behavior of particle systems 
in random environments seems to be a natural step ahead. A simple example 
of a single-file diffusion with disorder is what we call the exclusion process with 
variable diffusion coefficient. In this model, to each particle we attach a different 
diffusion coefhcient, and we let the system evolve in such a way that each particle 
performs a continuous- time, simple random walk with its own diffusion coefficient, 
conditioned to have at most one particle per site (the exclusion rule) . For a physical 
interpretation of this model and applications, see [3], [1]. 

In these notes we generalize the results of [2] , [5] to the simple exclusion process 
with variable diffusion coefficient. We will see that the scaling limits depend on the 
diffusion coefficient only through an effective diffusion coefficient, obtained from 
the related homogenization problem. Due to the variable diffusion coefficient, the 
approach of [2], [5] can not be repeated here, since in the associated stirring process, 
particles not only exchange their positions, but also their diffusion coefficients. We 
adopt here the approach of [IB], [H]. There, scaling hmits for the tagged particle 
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are related to scaling limits for the current through the (—1,0) bond in a zero-range 
dynamics. We call this relation the coupling method. The main problem is that 
the associated zero-range dynamics in the case of variable diffusion coefficient is a 
particle system with bond disorder. As proposed by Rost and Vares |19] . scaling 
limit for the current can be obtained from the hydrodynamic limit and density 
fluctuations of the density of particles. To obtain those scaling limits in the presence 
of random media, various methods have been developed [7], |18j . [6]. However, the 
method developed recently in [5] applies directly to our situation, avoiding the 
hard non-gradient tools needed in the other references. This is specially convenient 
in this setting, since the so-called spectral gap for the zero-range process in not 
uniform in density for our model. In fact, we think of this work as a nice direct 
application of the results in [8J. 

This article is organized as follows. In section [2] we define the model and state 
the main results. In section[3]we construct the associated zero-range dynamics and 
we prove the main results. In section |4] we discuss the relevance of the assumptions 
made on initial conditions and on the variable diffusion coefficients. 

2. Definitions and results 

Consider a system of continuous-time, interacting random walks {xi{t); i G Z} on 
the integer lattice Z, and a sequence of positive numbers {A^; i G Z}. The dynamics 
of this system is the following. The particle Xi {t) waits an exponential time of rate 
Ai, at the end of which it chooses one of its two neighbors with equal probability. 
If the chosen site is empty, the particle jumps to that site; if the site is occupied for 
another particle, the jump is not realized and a new exponential time starts afresh. 
We call this interaction the exclusion rule. This dynamics corresponds to a Markov 
process defined on the state space fig^ = {x G Z^; 7^ for i ^ j} and generated 
by the operator 

L§V(x) - ^'{l^^ + ^« ^ ^0") [/(x + e.) - /(x)] 

+l(x-e,Gl}g^)[/(x-e,)-/(x)]}. 

The number Ai is interpreted as the diffusion coefficient of particle Xi. We 
call the system {xi{t)}i described in this way, the simple exclusion process with 
variable diffusion coefficient, and we call {Ai}i the diffusion coefficient of {xi{ty}. 
Notice that when Ai = A does not depend on i, the process {xi{t)} corresponds 
to a labeling of the usual simple exclusion process. We refer to this situation as 
the homogeneous case. It will be of interest to define the particle configuration 
r]t G Vl'^^ = {0, 1}^ by r]t{x) = 1 if Xi{t) — x for some i ^ "L. Notice that 77* is not 
a Markov process, since in this case particles are indistinguishable, and therefore 
we can not keep track of the different diffusion coefficients. It will be useful as 
well to define the environment process r]f™{x) = r]t{x + xo{t)). Differently from 
the process rjt, 77!"" is indeed a Markov process. In fact, due to the exclusion rule 
and the nearest-neighbor jumps, the relative ordering of particles is preserved by 
the dynamics. From now on, we assume, relabeling the particles if needed, that 
Xi{t) < Xi+i{t) for any i G Z. Therefore, in the process r]'^"'"{t), the i-th particle at 
the right of the origin has diffusion coefficient A.; , and analogously for particles to 
the left of the origin. 
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Our aim is to study the asymptotic behavior of a tagged particle for this system. 
Let us assume that a;o(0) = 0, that is, that the 0-th particle starts at the origin. 
In the homogeneous case, it is well known that the behavior of a tagged particle is 
sub-diffusive, with a scaling exponent equal to 1/4 (instead of the diffusive exponent 
1/2). We will sec that this is also the case for the exclusion process with variable 
diffusion coefficient. For simplicity, we will assume that the sequence {Xi}i is an 
i.i.d sequence of random variables satisfying the ellipticity condition eo < Xi < Cq 
for any i S Z. More general conditions on {Ai}i will be discussed below. 

Now we explain how do we construct the initial configuration for the process 
{xi{t)}. Let Wo : R ^ [0, 1] be a continuous function. For each n > 0, we define the 
initial distribution /i„ as the product measure in fi*^^ with marginals 

Mn{??(a;) = 1} = 

We define a;o(0) = 0. For i > 0, we define Xi(0) recursively by Xi{0) = inf{x > 
Xi-i{0); r}{x) = 1}. For i < 0, we adopt a similar procedure. We denote this random 
initial configuration by -x.^^. Now we are ready to state our first result, which is a 
law of large numbers for xq {t) . For technical reasons, we will restrict ourselves to 
the case on which hma;_^±oo uo{x) = uo(±oo) both exist and are different from 0, 1 
and on which sup^, ua{x) < 1. 

Theorem 2.1. Let {xi{t)} be the exclusion process with variable diffusion coeffi- 
cient {Xi}i and starting from x"^. For almost every realization of the environment 
{Xi}i, 

n— *oo n 

in probability, where r{t) is the solution of 

r{t) .t 

u{t,x)dx = — / dxu{s,0)ds, 
Jo 

u{t, x) is the solution of the hydrodynamic equation dtu = Xd^u with initial condi- 
tion uq and X = {E[X~^])~^ . 

Notice in the previous theorem, the subdiffusive scaling, the dependence on 
the effective diffusion coefficient A and the quenched statement about the random 
diffusion coefficients. 

Once we have obtained a law of large numbers, a natural question is whether a 
central limit theorem holds. 

Theorem 2.2. Consider an initial flat profile uo{x) = tp. Let {xi[t)} be the exclu- 
sion process with variable diffusion coefficient {Xi}i and starting from x^. Then, 



1 2 / (l-y)va ^ 
hm —-prXo{tn ) = \ Z, 

n— >oo n / y 



■t, 



where Zt is a symmetric fractional Brownian motion of Hurst index 1/4 and the 
convergence is in the sense of finite distributions convergence. 

Notice in the previous theorem the explicit dependence of the limiting process 
in A and 
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3. The coupling method 

As we mentioned before, the environment process is a Markov process, but 
with a complicated non-local dynamics due to the ordering of particles. In order 
to obtain a local dynamics, we define a new process with state space il^^ = Ng 
by taking ^t(i) = Xi+i{t) — Xi{t). It is easy to verify that is a Markov process, 
generated by the operator 

L'^fio = E ^^{aim) [fic-'^') - /(?)] + + 1)) [fie^'n - /(o] }, 

where 5 : No ^ M is defined by g(0) = 0, g{i) = 1 if i > 1, and ^^'^ is defined by 

r'^(0 = u(j) + i, i^j 

The process is known in the literature as a zero-range process with bond dis- 
order. Although this functional transformation of the simple exclusion into a zero- 
range process was already exploited by Kipnis [14] to obtain the scaling limit of a 
tagged particle in the asymmetric simple exclusion process, the mathematical tools 
in order to treat the zero-range process with bond disorder have been developed 
only recently [8], [12]. The dynamics of this model can be described as follows. 
At each bond (i, i -I- 1) we attach a conductance Xi, corresponding to a exponential 
clock of rate A^. Each time this clock rings, we choose one of the directions i ^ i-{-l 
or i + 1 — > i with equal probability. Then we displace a particle from the origin 
site to the destination site. If there is not particle on the origin site at that time, 
nothing happens. We do this independently for any bond (z, i H- 1). 

Let us denote by jx.x+i{t) the total current of particles through the bond {x, x+1) 
for the process ^4. The key observation is that the position of the tagged particle 
XQ{t) in the original model corresponds to the current through the bond (—1,0) in 
the model Therefore, Theorem 12. II is consequence of the following result: 



Theorem 3.1. Under the assumptions of Theorem \2.1\ 
1 



lim o(i?i^) = - / \^{v{s,{)))dxv{s,Q)ds, 

n-*ca n ' JO 

where v{t,x) is the solution of the Cauchy problem 

f dtv{t, x) = dx (a$ {v{t, x))dxv{t, x)) 

[i;(0,x) =^0(2;), 



(3.1) 



the function ^{v) is given by $(w) = 1/(1 -|- w)^ and the initial condition vo{x) is 
defined below. 

About the initial condition vq, we notice that our choice of x^^ induces a choice 
in the initial distribution of as well. Therefore, we define Vq as the functional 
transformation of Uq in the continuous limit. This is given by the relation 

z+f^^ v{x)dx 

u{y)dy = z. 
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Defining V{z) = z + vo{x)dx, we see that V satisfies the ODE 



\V{0) 



0, 



from where we obtain vq by taking vo{x) — V {x) — 1 = (l — uo{V{x))) /uo{V{x)). 
In the same way, Theorem l2.2l is consequence of the following result: 



Theorem 3.2. Under the assumptions of Theorem 



1 2 

hm — jj_i,o(in ) = \ r Z, 

3.1. The law of large numbers. Equation 13. II corresponds to the hydrodynamic 
equation of the process ^t- For the process restricted to a finite volume, the 
so-called hydrodynamic limit has been obtained in [S] (see also [H]). Under the 
technical condition w(±oo) ^ 0,1, the result in [8] can be extended to infinite 
volume as in [4]. Therefore, we have the following 

Proposition 3.3. For any continuous function G : M ^ M of bounded support, we 
have 

lim —'S^^tn2{x)G{x/n)= [ G{x)v{t,x)dx 

in probability with respect to the law of S^t starting from , where v{t,x) is the 
solution of the hydrodynamic equation \3.1l 

Starting from this proposition, the law of large numbers for j-i,o(^) is easy to 
understand. In fact, the macroscopic density flux across in the hydrodynamic 
equation is given by Jo{t) ~ — \^{v{s,0))dxv{s,0)ds. Therefore, Theorem 13.11 
states that the current j_i,o(^), when properly rescaled, converges to Jo{t)- Notice 
that Proposition l3 . 31 does not give the desired result immediately, since an exchange 
of limits between the derivative and the scaling has to be justified. In order to 
justify this exchange, we can proceed as in |11| . so we will just give an outline of 
the arguments. When the total mass J vo{x)dx of the initial condition vo{x) is 
finite, the integrated current satisfies 

/>oo 

Jo(t) = / {v{t,x) ~v{Q,x)}dx. (3.2) 
Jo 

The identity is just saying that the total current across is equal to the mass at 
time t to the right of the origin, minus the mass at the right of the origin at time 0, 
an integrated version of the conservation of mass. Let us define Go(cc) = l{x > 0}, 
the Heaviside function. We can rewrite Jo{t) = J^Go{x)(^v{t, x) — v{0,x))ds and 
in an analogous way, 

-j-i,o{tn^) = - V{6„2(x)-Co(a;)}Go(x/n). (3.3) 

If the initial profile has non-null asymptotic densities t;o(±oo) = lima;^±oo vo{x), 
then formula (|3.2p still makes sense, but we can not send ri to oo in p.3p . since 
Go is not of compact support. Let us define Gi{x) = (1 — x/l)^ for a; > and 
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Gi{x) — for X < 0. Here (•)+ denotes positive part. Taking the Cauchy limit of 
the right side of the foUowing expression, it is not difficuh to see that 

-j-i.o{tn^) = lim - V {6n2(a;) -^o(a;)}G';(a;/0, 

the Hmit being in £^(P„), where P„ is the law of the process ^t„2 starting from . 
Moreover, the limit is uniform in I. But now, since the support of Gi is compact, we 
can use the result of Proposition [331 to pass to the limit in the previous expression, 
and to obtain the desired result. 

3.2. The central limit theorem. We recall the strategy in order to prove a law 
of large numbers for the current j_i_o(i) followed in the previous section. The 
idea is to relate the current to some linear functionals of the density of particles. 
This relation holds for the microscopic model, and also for the density flux in the 
macroscopic equation. In this way, the scaling limit of the current should follow 
from the scaling limit of the density of particles. However, the functional considered 
was not continuous, and therefore a cut-off argument was needed in order to take 
the limit. The final step was to justify an exchange of limits between the scaling and 
the cut-off. This strategy also gives a central limit theorem for the current in the 
zero-range process ^t- Differently from [llj . in the case of the zero-range process 
the central limit theorem for the density of particles has only been obtained in 
equilibrium, restricting ourselves to that situation. 

Now we explain what do we mean by a central limit theorem for the density of 
particles. For the initial configuration x^, the associated initial distribution is 
a product, uniform measure in with marginals given by 

A simple computation shows that the measure satisfies detailed balance with 
respect to the evolution of ^t, and therefore it is an invariant, reversible measure 
for this dynamics. It is not hard to check that it is also an ergodic measure for ^t. 
Let us notice that the number of particles in the zero-range process is locally 
conserved, and due to ergodicity, it is the only locally conserved quantity. Therefore, 
it is more natural to parametrize the family \v^\Lp S [0,1]} of ergodic, invariant 
measures by the density of particles at each site. It is easy to see that the density 
of particles p(}p) defined by p(}p) — f £,{x)D^ satisfies p{ip) = (1 — ip)/ip. We 
define then ip{p) = 1/(1 -I- p) and i>p = v^p^^py Now fix a density p € (0, cx)) and 
consider the process ^t, starting from the measure Vp. This is equivalent to start 
the exclusion process with variable diffusion coefficient from the configuration x" 
for V3= 1/(1 + p). 

Let us define the 5'(R)-valued process 3^" by 

(G) = ^ ^ G{xln) (^,„. (x) - p) , 

where G is in 5(]R), the Schwartz space of test functions in M. We call this process 
the density fluctuation field. The following proposition has been proved in [8] : 

Proposition 3.4. The fluctuation field converges in the sense of finite distribu- 
tions to a generalized Ornstein- Uhlenbeck process of mean zero and characteristics 
\^{p)dl and ^\p/{l + p)d^. 
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As in the case of the law of large numbers, the current j-i^o{t) satisfies 
^j_i,o(in2) = lim y^iGi)~ySiGi), 
the limit being uniform in £^(P„). Therefore, we conclude that 
lim J-j_^^o{tn^) = lim yt{Gi)-yoiGi), 

n^oc l—i-oo 

from where Theorem 13.21 follows. 

4. Conclusions and remarks 

1. Comparing our results with the classical result of Arratia [2l and the results in 
[T2] , we see that the variable diffusion coefficient enters into the scaling limit of the 
tagged particle only through the effective diffusion coefficient A. Notice that this 
effective diffusion coefficient is the same obtained by homogenization, although our 
system is subdiffusive. This is due to the fact that, for particle systems with bond 
disorder, the disorder and the time evolution can be decoupled using an stochastic 
version of the compensated compactness lemma [S] . We conclude that the subdiffu- 
sive behavior is a purely dynamical feature, since the factor \J (\ — (f)\^/(p is the 
same appearing for a simple exclusion process with uniform diffusion coefficient A. 
In particular, the variance of XQ{t) is of order ^y2Jd/^^{^ — ip) j ip a& t ^ oo. 

2. About our assumptions on the initial conditions for the law of large numbers, the 
product structure of the initial measure /i„ is far from necessary. The more general 
assumptions are the following. First we need the initial measure to be associated 
to some macroscopic profile uq, in the sense that the particle density should converge 
in distribution to uo(-) (for a precise definition, see Section 2.3 of [8]). The second 
assumption is a bound on the entropy density of the initial measure /z„ with respect 
to the measure /i*^ obtained as the product measure associated to the flat profile 
uq = if. We also need the distribution of to be stochastically dominated by the 
measure Vp for some p > 0. This imposes the restriction sup 2, 1*0(2^) < 1- 

3. In the CLT for the tagged particle, a key input is the fluctuation result for the 
particle density. Proposition 13.41 A complete proof of the results in 13.41 is only 
available in equilibrium, making the choice of yH^ rather restrictive. To extend 
Proposition 13 . 41 to the non-equilibrium setting remains one of the main open prob- 
lems in the theory of scaling limits of particle systems. Notice that such a result 
have been obtained in |12j for the simple exclusion process in random media and 
in [To] for the zero-range process. None of these approaches can be adapted to our 
situation, however. 

4. The assumptions on the variable diffusion coefficients are not, by far, 
optimal. We have taken as an i.i.d. sequence just to fix ideas. In fact, in 
order to get the results of Theorem 12.11 it is enough that sup, A^ > and 

1 ° 1 " 

lim - y Ari = lim - V A^' = A-^ 

i— — n i—0 

both exist and are different from 0. To obtain the results in Theorem 12. 2[ we need 
in addition the ellipticity condition e^^ > Ai > cq for every i. In particular, both 
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results remain true if for instance we take a periodic sequence {Xi}i of diffusion 
coefficients. 
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